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Stability and instability bands in classical mechanics are well-studied in connection with systems 
such as described by the Mathieu equation. We examine whether such band structure can arise in 
classical field theory in the context of an embedded kink in 1+1 dimensions. The static embedded 
kink is unstable to perturbations but we show that if the kink is dynamic it can exhibit stability in 
certain parameter bands. Our results are relevant for estimating the lifetimes of various embedded 
defects and, in particular, loops of electroweak Z-string. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



A wide variety of classical static solutions in field theo- 
ries have been constructed and they can play crucial roles 
in particle physics, condensed matter systems, and cos- 
mology. Commonly the static solutions exist because of 
non-trivial topology of the field theory and in these cases 
the solution is known as a "topological defect" . Alter- 
nately the solution may exist in the form of a topological 
defect embedded in a theory where the requisite topology 
is absent. In these cases the static embedded topological 
defect is usually unstable and can decay. We are inter- 
ested in determining if this instability is also present in 
embedded defects that are not static. 

It is helpful to consider a mechanical analog of the 
unstable defect solution. Consider a particle in two spa- 
tial dimensions with position x(t) = (x(t),y(t)) and La- 
grangian 



where 



1 
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with k x , k y and c being positive parameters. The ex- 
ternum of the potential is at x = and is a saddle point. 
The positive eigenvalue at the saddle point is along the 
x direction and the negative eigenvalue is along the y 
direction. The equations of motion are: 



mx 



c x — cxy 



my = +k y y 



cx 2 y 
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where overdots denote derivatives with respect to t. The 
static solution is x = = y and is clearly unstable to 
perturbations in the y direction. This is seen by taking 
y to be a perturbation and then linearizing the equation 
for y: 



This equation has exponentially growing solutions and 
hence the particle is unstable to "rolling off" in the y- 
direction. 

Next consider a dynamic solution of the system: 



o(t) = Acos(ujt) , y(t) = 



(6) 



where uj = ^k x /m. Is this solution stable? 

Once again we proceed by perturbing the y equation, 
and linearizing in the perturbations. This gives 



my = +k y y — cA 2 cos 2 (ujt)y 



(7) 



where y denotes the perturbation. This is the Mathieu 
equation (Sec. 5.2, 0). 

There exist solutions to the Mathieu equation of the 
form 



y(t) 



(8) 



a result which is derived using Floquet's theorem. The 
value of v depends on the parameters k y /k x and cA 2 jk x 
and can be found as described in Sec. 5.2 of Ref. Q. If 
v is imaginary, the original solution is unstable; if v is 
real, the original solution is stable. 

A remarkable feature of this system is that there are 
bands of stability and instability. For fixed w, if A is 
close to zero, there is an instability and the particle rolls 
off in the y-direction. For larger A, we find stability - 
the particle gets displaced in the y-direction at the saddle 
point, but is unable to roll off by the time it has moved 
away from the saddle point in the x-direction. This is 
shown in Fig. |] where we have plotted the stability and 
instability bands in the a — q plane, where 
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Then, for fixed k y /k x , we get the line 
a=-2 q - k f 



(9) 
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on the plot. This line intersects the unstable and stable 
regions alternately as we go to larger amplitudes for fixed 
model parameters. 
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FIG. 1: Band structure for the Mathieu equation. For q < 0, 
the stability bands lie between consecutive dashed or consecu- 
tive solid curves. The straight line denotes a set of parameters 
for the mechanical model for different values of the oscillation 
amplitude. Larger amplitudes correspond to points on the 
line that are further away from the origin. As the amplitude 
is increased, the line increasingly lies in the stability bands but 
no matter how large the amplitude, there are always bands 
where there is instability. 



It is tempting to think that the pattern of stability 
and instability as the amplitude is increased can be un- 
derstood in terms of the time spent by the particle in 
the vicinity of the saddle point as compared to the de- 
cay time of the static solution. However, this intuition 
does not work since for yet larger A we again get insta- 
bility. The band structure then has to be understood as 
a resonance phenomenon between the various oscillatory 
modes of the particle. 

If we include frictional forces on the particle as it moves 
on the potential, the amplitude of oscillations will gradu- 
ally diminish and eventually the particle will enter a band 
of instability. The decay time of the oscillations will then 
be determined by the time it takes for the particle to slip 
from a stability band to an instability band. 

We are interested in determining if a similar phe- 
nomenon can occur in classical field theories. If a cer- 
tain static field solution is unstable, can its lifetime be 
much longer due to some dynamics? There is clearly an 
affirmative though trivial answer - since these are rela- 
tivistic field theories, a boosted field solution lives longer 
due to relativistic time dilation. However, we are not 
interested in this factor and would like to determine if 
there is an effect similar to the band structure observed 
in the mechanical problem. There is another complica- 
tion in field theory: the dynamical solution will almost 
certainly emit radiation and dissipate and hence is simi- 
lar to the mechanical case with the inclusion of friction. 
Yet in the examples we will consider, the dissipation due 
to radiation is very small and can be ignored. 

We will restrict our attention in this paper to classical 
field theory in 1+1 dimensions. Our first approach to the 
problem is to construct unstable (embedded) kink solu- 
tions and study the growth of linearized perturbations on 



kink solutions that are forced to oscillate. The similar- 
ity between this field theory problem and the mechanical 
problem described above can be understood by going to 
the rest frame of the kink as it oscillates. In this (non- 
inertial) frame the perturbation modes are oscillating in 
the vicinity of the kink. In the time when a perturbation 
mode lies outside the kink, it cannot grow. The mode 
can only grow during the time that it lies within the 
kink. This is exactly like the mechanical model where 
the particle (analogous to the perturbation mode) can 
roll off in the y-direction near the saddle point but not 
when its oscillations take it away from x = 0. Hence 
we expect the behaviour of the perturbation mode to be 
similar to that of the particle. 

Indeed our analysis reveals band structure in the 
growth of the perturbation modes though we only find 
stability when the amplitude of the kink oscillations is 
comparable to the width of the kink. In other words, in 
the linearized analysis only "jittering" kink backgrounds 
are found to have stability bands. We have also per- 
formed a second analysis of the problem. In this ap- 
proach, we impose an external potential in which the 
kink can oscillate. Now we evolve the full field equations 
and find the lifetime of the kink. Here again we find band 
structure. Furthermore, the amplitude of motion can be 
significantly larger than the kink width. We attribute 
the enhanced stability bands at large amplitudes to the 
reduced and varying kink width in the external potential. 

Our study has obvious implications for the study of 
classical solutions that are unstable when static. Ex- 
amples of such solutions include a variety of embedded 
defects ^ ||. In particular, electroweak Z-strings are an 
important example |^|, |5|. Assuming, for now, that the 
band structure exists for dynamical solutions in 3+1 di- 
mensions as well, we discuss the possible relevance of our 



work to the electroweak model in Sec. IV 



II. 



LINEARIZED ANALYSIS OF FIELD 
THEORY PROBLEM 



The 1+1 dimensional field theory we choose to study 
is one that contains embedded kink solutions ^, || . The 
complex scalar field, <f> — <p\ + i<f>2, has the Lagrangian 
density: 



L=\\d ix <t>\ 2 -V{<j>) 



(11) 



where 



v^) = \m 2 -i? + K^m- as) 

(In 1+1 dimensions, the fields arc dimcnsionless. For con- 
venience, we have also rescaled the coordinates so that 
they are dimensionless.) The last term in the potential is 
analogous to the cross-term in eq. (J2j) and it violates the 
U(l) symmetry <fi — > exp(ia)<fi when c ^ 1 but preserves a 
Z2 x Z2 subgroup. A vacuum expectation value of either 
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the real (0i) or imaginary component (02) of spon- 
taneously breaks the symmetry to the identity group. 
Hence there are topological kink solutions in the model 
for c 1. However, we will not be interested in these 
topological solutions. Instead we will study the embed- 
ded kink solution which is not topological and exists for 
all values of c. This solution is: 



where 



( x 

f>i = tanh — = 
VV2 



= 



(13) 



In the sub-space of field configurations defined by 02 — 
0, the embedded kink is topological and hence stable. 
However, in the full field space, the static embedded kink 
is unstable to perturbations in the 02 direction. The 
equation for linearized perturbations in the 02 direction 
in the static embedded kink background is: 



- 02 + (O 



D02 = 



(14) 



where 0i is given in eq. (|13|). The solution for 02 is a 
hypergeometric function (Sec. 6.3, E|): 



h = e-sech" (^) 

f(k + ^ + P,K+^-P\K+1\z) (J.r,) 



with 



1 

Z ~ 2 



K 



1 + tanh I -^= ) 
V2J 



2(c- 1 



P = w2c+i . 
The lowest eigenvalue vq is given by: 



1 L 1 5 
2V 2c+ 4-4 



(16) 



(17) 



(18) 



(19) 



Since Vq is imaginary for 1 < c < 3, the static solution 
in eq. ( p"3] ) is unstable for these parameter values. For 
c = 3, the lowest value of v is zero, corresponding to the 
translation mode of the topological Zi kink (fj). Hence 
for c > 3, the static embedded kink is stable. 

We now want to introduce some dynamics to the em- 
bedded kink solution and study the effect of the dynam- 
ics on the unstable mode in the 02 direction. In the 
linearized approach to the problem, the dynamics is in- 
troduced "by hand" by introducing a time dependence 
in the unperturbed background: 



b\ = tanh 



= 



(20) 



g(t) = A cos(wi) 



(21) 



describes oscillations with amplitude A and angular fre- 
quency w, and 7 is the "Lorentz factor" 



1 



7 



r 



(22) 



We require ojA < 1 so that the speed of the kink is al- 
ways sub-luminal. In this background the equation for 
02 perturbations is: 



(c tanh 2 A - 1)02 = (23) 



where 



X(t,x) = -L(x-g(t)) . (24) 

Note that boosts of the static solution can be obtained 
by setting g — vt where v is the speed of the kink. Such 
backgrounds do have a longer lifetime but this is due to 
the special relativistic time dilation. 

In Appendix A we extend Floquet's theorem to field 
theory to show that solutions of eq. (23) have the form: 



<j> 2 {t,x) = e lut P{t,x) 
where P{t 1 x) is a periodic function in t: 



(25) 



P(t + 2hlj- 1 ,x) = P(t,x) 

The next step is to find the eigenvalue v. We were unable 
to determine v analytically and instead solved eq. ( |2^ ) 
numerically for various parameters A, uj and c and with 
an initial perturbation: 



02 (t = 0,x) = sech(X(0,a;)) 



(26) 



A sample result of the numerical evolution is shown in 
Fig. ^. Here we plot P(t, x) for several different values 
of X (see eq. (p4|) ) versus time. The plots are consis- 
tent with the form in cq. (p5|). For some values of A, 
02 (i, x) grows exponentially with time, while for others 
it oscillates without growth. This allows us to find the 
lifetime, r, of the background. In Fig. ||[ we show the 
inverse lifetime, r , as a function of the amplitude A for 
c = 1, 2, 2.5 with ui fixed by the relation u = 0.99/ A. The 
reason we held uj = 0.99/A is that we expect that sta- 
bility will depend on the acceleration a c of the kink and 
this is < uj 2 A. Since the speed cannot exceed the speed 
of light, we impose luA < 1. Therefore a c < v 2 lax /A 
where v rnax = luA is the maximum speed of the kink. So 
we expect greatest stability when v ma x ~ 1 and for small 
oscillation amplitude. In Fig. the bands of stability 
occur where the inverse lifetime vanishes. 

It is worth noting that the values of A for which sta- 
bility can occur are of order the half-width of the kink 
which is given by Hence, stability bands only oc- 
cur when the background is undergoing a rigid jittering 
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III. FIELD THEORETIC PROBLEM WITH 
EXTERNAL POTENTIAL 
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FIG. 2: The spatial value of the function P(t,x) for X = 
(solid curve), 1 (dashed curve), 2 (dotted curve) as a function 
of time. 
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FIG. 3: The decay rate r _1 is plotted versus the amplitude 
of oscillation for c = 1 (solid curve), 2 (dashed curve), and 
2.5 (dotted curve), as determined in the linearized analysis. 
Stability bands are regions where r _1 = 0. There are no 
stability bands for c = 1 but these do exist for c = 2 and 2.5. 
The curve for c — 2 has been off-set by 0.15 along the y-axis 
for clarity. 



motion. The point A — is a singular point in our ap- 
proach since the oscillation frequency is to = 0.99 /A and 
this becomes infinite as A tends to zero. Another im- 
portant point is that there are no stability bands for the 
embedded kink in the U(X) model obtained with c = 1. 
We also observe that the dynamic kink with c > 3 can 
show instability. This is not inconsistent with our earlier 
analysis of the static kink - only the static kink is sta- 
ble for c > 3 and there is no reason why rigidly moving 
dynamic kinks should be stable when c > 3. 



A drawback of the linearized analysis is that it keeps 
the background "rigid" - the profile function is always 
a hyperbolic tangent - while an oscillating kink should 
also experience some profile fluctuations and be more like 
a "fluid" . Hence in our second analysis of the problem 
we treat the full field theory, rather than perform a per- 
turbative analysis about a chosen field background. To 
obtain a dynamic embedded kink, we have introduced a 
background potential in which the embedded kink can 
oscillate. The Lagrangian density now is given by eq. 
( pi] ) together with an additional "external" term in the 
potential: 



v ext {<t>) = \u{x)^\-iy 



(27) 



where U (x) is a potential well in which the embedded 
kink oscillates. The external potential can possibly orig- 
inate due to the expectation value of other very massive 
fields on which backreaction is unimportant 1 . A sim- 
ple example for the external potential is U(x) = dx 2 /2 
(d > 0). Then V ext provides a force between the posi- 
tion of the embedded kink (inside which fa ~ 0) and 
the origin of the coordinate system (x = 0). Hence an 
embedded kink that is initially displaced from x = 0, os- 
cillates about x = and the frequency of oscillation is 
determined by the parameter d. In our numerical analy- 
sis, we have chosen 



U(x) = -a 2 (l 



-(a/a)' 



) 



(28) 



and we have chosen a = 15, d = 0.1. Then U(x) has 
the x 2 form for small x and changes to a constant for 
x >> a. The feature that U(x) does not diverge at large 
x has the advantage that radiation of fa waves from the 
oscillating kink can escape to infinity. 

Now the equations of motion for fa and fa are: 

fa - fa{ + (fa 2 + aj>\ - I) fa + 2U(fa 2 - I) fa = (29) 



fa - fai + {cfa 2 + fa 2 



i) 







(30) 



The equations of motion are first solved numerically 
with the initial condition: 

fa(t = Q,x) =tanh(X(0,a;)) , fa(t = 0,x) = Q. (31) 



1 This could be realized physically in a system which has a phase 
containing only unstable kinks (call it the A phase) that is sand- 
wiched between layers of another phase (call it the B phase) in 
which the kinks are stable. Then the B-phase stable kink could 
be driven to oscillate across a region of A-phase where it would 
be unstable. The potential V e xt contains both the driving poten- 
tial (U(x)) and also the spatial variation of the phases because 
of the field dependence in it. 
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IV. CONCLUSIONS AND DISCUSSION 
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Amplitude 



FIG. 4: The decay constant r _1 as a function of the initial 
amplitude of oscillation in the full field theory for c = 1 (solid 
line), 2 (dashed line), 2.5 (dotted line). Here too we observe 
bands of stability though none occur for c = 1. 



where 



X(0,x) = y/l + 2U(A)X(0,: 



(32) 



The factor of 1 + 2U(A) arises since, as seen in eq. (|29|) , 
the external potential contributes to the width of the 
kink. On evolution the amplitude of oscillations is ex- 
pected to decay due to radiation. However, the decay is 
very slow and we can easily track the evolution of the 
fields through several tens of oscillations with the am- 
plitude of oscillation nearly constant. Another feature 
worth pointing out is that the angular frequency of os- 
cillations depends on the amplitude of oscillation. The 
angular frequencies of the runs with large initial ampli- 
tudes are ~ 1/A but at small amplitudes the angular 
frequency saturates at ~ 0.3. 

Next we ran our numerical code with a non-vanishing 
perturbation. The initial condition for <pi is still given 
by eq. ([u]) but 2 is chosen to be given by 



02 (t = 0,x) = esech(AT(0,x)) 



(33) 



where e is a very small number. The numerical analysis 
is done with fixed external potential. 

As in the linearized analysis, a plot of <f>2(t,x) versus 
time either shows exponential growth (superposed on os- 
cillations), or only oscillations. This allows us to calcu- 
late the rate of growth of the perturbations. Once again 
we find bands of stability (see Fig. |J). The new feature is 
that these bands occur even for rather large amplitudes 
showing that dynamical (and not only "jittering" ) em- 
bedded kinks can be stable. Note that stability bands 
do not occur in the U(l) model (c = 1) just as in the 
linearized analysis. 



Our analysis of the dynamical embedded kink has 
shown that there are stability bands in field theories, sim- 
ilar to those for the Mathieu equation. This means that 
for certain dynamics the embedded kink is stable to per- 
turbations. In these cases, the decay of the embedded 
kink is determined by the rate of radiation and hence its 
lifetime can be much longer than the lifetime determined 
by considering the decay of a static embedded kink. 

It is worth commenting on the differences between the 
linearized and full field theory approaches since the re- 
sults are quite different. In the linearized case, the os- 
cillating background is inserted by hand and this gives 
stability bands occurring when the kink has the largest 
acceleration. In the full field theory, an external poten- 
tial is imposed, supposedly arising from an external field. 
Now the kink oscillates naturally, but we observe stabil- 
ity bands at large amplitudes. This feature is explained 
by our earlier remark that the external potential influ- 
ences the structure of the kink. As the kink gets further 
away from the center of attraction, it gets thinner and 
hence more stable. (The connection between the width 
of a static kink and the instability eigenvalue can be de- 
rived quite easily from a generalization of the solution in 
eq. (|5|) as given in Ref. M.) Then the reason why kinks 
with large oscillation amplitudes show more stability is 
that they spend more time in regions where the kink is 
thin and relatively stable. This is in addition to the ef- 
fect described in the introduction: the unstable mode 
moves in and out of the kink, and can grow only while 
it lies within the kink. Yet the novel feature of stability 
and instability bands must still be understood in terms 
of resonance between the oscillation frequency and the 
phase of the unstable mode when it enters the kink. 

A natural question at this point is - can the tension 
forces that cause extended defects such as domain walls 
and strings to oscillate be described by an external po- 
tential of the type we have used? The answer is in the 
negative, since the tension forces do not cause changes 
in the thickness of the defects. The only reason for the 
defect thickness to change is due to Lorentz contraction. 
This is our reason for believing that an extended defect 
will be more accurately described by the linearized anal- 
ysis of the previous section and not by the full field the- 
ory analysis. We are currently extending our analysis 
to higher dimensions with the aim of studying extended 
defects. 

In the particle physics context, the results are most 
relevant for the Z-string solutions in the standard model 
of the electroweak interactions. In 1977 Nambu @ had 
estimated the lifetime of segments of Z-string ("dumb- 
ells") by calculating the electromagnetic emission from 
the magnetic monopoles at the ends of a rotating dumb- 
ell. If, however, we have a closed loop of Z-string, 
there are no magnetic monopoles that can lose energy 
by emitting electromagnetic radiation and the only de- 
cay processes are by radiation of massive Higgs and gauge 
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bosons, and by the structural instabilities of the Z-string 
discussed in Ref. ]§| . (These may be understood in terms 
of W-condensation || in a Z magnetic field background 
ju| 0, flf, pi.) If we assume that the Z-string in 3+1 
dimensions also exhibits stability bands, it would mean 
that the lifetime of certain Z-string loops may be deter- 
mined by their radiation rate Jl4| and not by the insta- 
bility growth rate of a static Z-string. This would yield a 
much longer lifetime and would open the possibility that 
Z-string loops might be produced in future accelerator 
experiments and might survive long enough to enable 
detection. 

Some cautionary remarks in extrapolating our results 
to Z-string loops are in order. The Z-string loop is an 
extended object and if it is to have an extended lifetime, 
every point on the loop should lie in a stability band. 
It may even turn out that the Z-string is like the c = 1 
embedded kink and that there are no stability bands for 
any dynamics. Furthermore, the dynamics of the string 
is due to tension and not due to an external potential. 
Hence the Z string stability problem could be like the lin- 
earized example, in which case only "jittering" Z strings 
might be stable. The Z string loop can also have angu- 
lar momentum and this will be an important factor in 
its evolution. (The stability of spinning string solutions 
has been studied in Em.) Hence we feel that it is impor- 
tant to directly investigate the lifetime of Z-string loops, 
keeping in mind the possibility of stability bands. 
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APPENDIX A 

The following derivation was motivated by Ince's 
derivation of Floquet's Theorem |lf|. Consider a dif- 



ferential equation of the form: 

4> + F(t,x)ct> = H<j) = 0. (Al) 

If H is Hcrmitian over an N dimensional Hilbert space, 
we can expand this as: 

N 

ff^ = 5^Ai0j(t,ar)=O (A2) 
»=i 

where the <j>i span the Hilbert space. Now, if F(t, x) = 
F(t + T,x) where T = 2tt/uj in the notation of the earlier 
sections, then for each solution (f>i(t, x), there is a solution 
(j>i(t + T, x) by the symmetry in H. Since the <j>i are 
complete, 

N 

4>i (t + T, x) = J] ajj (j>j (t,x). (A3) 
i=i 

Denoting the matrix [a^] by A and the N x N identity 
matrix by In, for each non-degenerate root Si of the 
characteristic equation 

det(A-sl N ) =0 (A4) 

we can construct a basis ipi from the 4n such that 

ii i (t + T,x)=s i il)i{t,x). (A5) 

Multiplying both sides of the above result by e~ tVi ^ t+T \ 
we see 

e -fi(*+rty i(t + Tj x) = s . e -iu iT [ e -<^.( t) x) ] . (A6) 

Then it is clear that e~ wt ipi(t,x) is periodic in t with 
period T if s,; = e +WiT . Therefore, solutions of the form 

rl> i (t,x) = e iu < t P(t,x) (A7) 

exist to the equation Hif>i = with P(t + T,x) = P(t, x) 
and vi given by e +-lViT = Sj for non-degenerate S{. 
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